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1 Introduction 



Abstract We show that in complete agreement with clas- 
sical mechanics, the dynamics of any quantum mechanical 
wave packet in a linear gravitational potential involves the 
gravitational and the inertial mass only as their ratio. In 
contrast, the spatial modulation of the corresponding energy 
wave function is determined by the third root of the product 
of the two masses. Moreover, the discrete energy spectrum of 
a particle constrained in its motion by a linear gravitational 
potential and an infinitely steep wall depends on the iner- 
tial as well as the gravitational mass with different fractional 
powers. This feature might open a new avenue in quantum 
tests of the universality of free fall. 



The equivalence principle is a cornerstone in the foundations 
of general relativity [1|. Indeed, the assumption of the pro- 
portionality of inertial and gravitational mass implies that in 
a linear gravitational potential all bodies experience the same 
acceleration and fall with the same rate. Without this univer- 
sality of free fall, the geometrization of gravitation and its 
reinterpretation as curvature of spacetime would not be pos- 
sible. The fact that several alternative gravitational theories 
predict the breakdown of the universality of free fall [2] is one 
of the main reasons that drives physicists to test this bedrock 
of modern physics to higher and higher accuracy. 

Motivated by the seminal papers on neutron interfero- 
metry [3-7 1 and the more recent, impressive matter wave ex- 
periments 1 8-14], we address in the present paper the ques- 
tion how the inertial and gravitational mass enter in non- 
relativistic quantum mechanics. We show that in the case of 



a linear gravitational potential, quantum dynamics only in- 
volves the ratio of the two masses in complete accordance 
with classical Newtonian mechanics. However, depending on 
the specific preparation of the initial state, inertial and grav- 
itational mass may appear in a more complicated way in the 
time evolution of a physical state. As an example of such an 
initial state, we discuss the energy eigenfunctions in a linear 
potential 1151161 . which have been analyzed e. g. in the con- 
text of the coherence of an atom laser IfTTl or in connection 
with the so-called atom trampoline, also known as the quan- 
tum bouncer [18-22]. Indeed, the energy eigenstate in this 
system is non-classical since the corresponding phase space 
equations for the Wigner function do involve Planck's con- 
stant h. As a result such states are ideal objects to study the 
role of inertial and gravitational mass in quantum mechanics. 

Three central results obtained in this paper stand out: (i) 
The quantum dynamics reduces to classical dynamics and 
therefore can only involve the ratio of the inertial mass rrii 
and the (passive) gravitational mass m g , (ii) the spatial mod- 
ulation of the energy eigenfunctions depends on the third root 
of the product of the two masses, and (iii) the energy eigen- 
values of the gravitational atom trampoline are proportional 
to (rrig/mi) 1 ^ 3 . 



1.1 Tests of the universality of free fall 

The universality of free fall, often referred to as "weak equiv- 
alence principle", states that all bodies experience the same 
gravitational acceleration independent of their internal struc- 
ture and composition, provided they are so small in size that 
one can neglect the effects of gravity gradients. In other 
words, the (inertial) mass of a body is proportional to its 
weight, with an universal proportionality constant. A viola- 
tion of this principle would arise e. g. when the interaction 
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energy between the nucleons in an atom would not contribute 
in the same manner to the gravitational mass, as it would for 
the inertial. 

The classical tests II 1121 of the universality of free fall as- 
sume that the specific gravitational acceleration g(A) does 
depend on the internal structure or the composition of the 
body A. This assumption translates in terms of the inertial 
and gravitational mass into the relation 



9(A)^g[^ 

mi 



(D 



Here, the gravitational acceleration g should be considered 
as a standardized acceleration corresponding to a particular 
reference body. 

A measure for the breakdown of the universality of free 
fall is the so-called Eotvos parameter 



T)(A,B) 



m a \ , ( m 9 \ 
mi j A \ mt J £ 



g(A) - g{B) 
g(A)+g(B)' 



which quantifies the normalized difference in the gravita- 
tional accelerations between two different bodies A and B. 

The first tests of the equivalence of inertial and grav- 
itational mass relied on pendulum experiments and can be 
traced back to Newton and Bessel [23|. A great step towards 
higher accuracies was realized by the classical torsion bal- 
ance experiments of Eotvos 11241 and Roll et al. (25). Cur- 
rently the best upper limits for the Eotvos parameter ll26ll 
come from lunar laser ranging on the one hand and from the 
so-called "Eot-Wash" experiment [27 28] on the other. The 
latter uses a sophisticated rotating torsion balance and limits 
the Eotvos parameter to 



ry(Be,Ti) = (0.3 ± 1.8) x 10 



-13 



for the gravitational acceleration of Beryllium and Titanium 
towards Earth. 

The motivation for quantum mechanical tests B29I30II of 
the universality of free fall stems from the increase in accu- 
racy that atom interferometry is expected to offer in the fu- 
ture. Matter wave interferometry with freely falling Rb 8j and 
Rb 87 isotopes has already been performed lfl2l . and several 
other experiments worldwide using different species of atoms 
are right now in preparation ll3"D . 



1.2 Discussion of related work 

In this paper we study the quantum mechanics of a particle in 
a linear potential. Needless to say, this topic appears promi- 
nently in many papers, in particular in connection with the 
atom trampoline |[T8lj20l whose energy eigenstates have been 
theoretically investigated in ETI . as well as in the context of 
cold neutrons [32]. 

The work closest to ours is that of Davies ll33l who has 
investigated the problem of a quantum mechanical particle 



in a linear gravitational field to gain insight into the equiva- 
lence principle. However, motivated by the classical motion 
his main emphasis is on a thorough analysis of the travel time 
of wave packets. In complete agreement with our conclusion 
that the dynamics is classical, he finds that classical and travel 
times agree far from the classical turning point. However, 
there are quantum corrections near the turning point. His cal- 
culations are based on the definition of the Peres quantum 
clock H341. 



1.3 Outline of the paper 

We start in Sect. |2]by recalling the universality of free fall in 
Newtonian mechanics. Since we are interested in the Wigner 
phase space formulation of the corresponding quantum me- 
chanical version, we first introduce in Sect. [3] the quantum 
Liouville equation describing the dynamics of the Wigner 
function in an arbitrary potential. Moreover, we present the 
partial differential equations in phase space determining an 
energy eigenstate in this potential. In Sect. [4] we then apply 
these equations to analyze the dynamics of the Wigner func- 
tion in a linear gravitational potential, as well as to determine 
the Wigner function of the corresponding energy eigenstates. 
Whereas the quantum dynamics just reflects the classical time 
evolution and does not depend Planck's constant H, the phase 
space analog of the energy eigenstates does display quantum 
features and involves K. This fact stands out most clearly in 
the energy eigenfunctions of the linear potential discussed in 
Sect. [5] where we show that the wave vector governing the 
spatial modulation of the probability density is determined 
by the third root of the product the inertial and gravitational 
mass. In addition, we examine the energy eigenfunctions and 
the eigenvalues of the atom trampoline, also known as quan- 
tum bouncer, that is a particle trapped in the bounded poten- 
tial resulting from the combination of a linear potential and an 
infinitely steep wall. We conclude in Sect. |6]by summarizing 
our results and by outlining possible experiments. 

In order to keep the paper self-contained we summarize 
concepts pertinent to the present discussion in several appen- 
dices. For example in Appendix [A] we recall that the time 
evolution of a particle in a linear potential can be represented 
in phase space as a product of a shearing operator followed 
by a displacement. It is only the displacement which contains 
the gravitational acceleration. This decomposition provides 
us with deeper insights into the physics of non-spreading 
Airy wave packets as outlined in Appendix iBl We dedicate 
Appendix [C] to a discussion of the semi-classical limit of 
the energy wave function in a linear gravitational potential. 
Within the Jeffrey s-Wentzel-Kramers-Brillouin ( J WKB) ap- 
proximation we regain the universality of free fall. Moreover, 
we can identify phase space quantization as the origin of the 
unusual scaling properties of the energy eigenvalues of the 
atom trampoline in terms of inertial and gravitational mass. 
We conclude in Appendix ITJ1 with a phase space analysis of 
the atomic fountain. 
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2 Universality of free fall in Newtonian mechanics 



2.2 Ensemble dynamics 



In the present section we briefly recall Newton's law of mo- 
tion for a single particle in a linear gravitational potential 
which shows that the dynamics only depends on the ratio be- 
tween the gravitational and the inertial mass. For this reason, 
the free fall of two particles of different compositions is iden- 
tical provided the inertial and gravitational mass are propor- 
tional to each other and the particles start with the same initial 
conditions. Note that the universality of free fall requires the 
mentioned proportionality constant to be independent of the 
particles composition, in which case it can be absorbed in the 
universal gravitational acceleration. 

Moreover, we examine the dynamics of a classical en- 
semble of identical particles falling in linear gravitational po- 
tential, thereby emphasizing some complications that come 
along with tests of the universality of free fall based on clas- 
sical statistical mechanics. 

Since we consider a homogeneous gravitational potential, 
we can restrict our analysis of the dynamics and kinemat- 
ics to one spatial coordinate z. Moreover, instead of dealing 
separately with the inertial and gravitational mass, we will 
take advantage of the particle dependent gravitational accel- 
eration (Q3 and denote the inertial mass of a particle most of 
the time by m = 



2.1 Single particle dynamics 

The time evolution of a particle moving in an external poten- 
tial V = V(z) follows from Newton's law of motion 



rrii z 



dV 

dz 



(3) 



where the dots indicate differentiation with respect to time. 
For a linear gravitational potential 



Vi(z) = m g gz 
we obtain from Eqs. (Q~|) and (0) 



z = 9 = -9- 

rrii 



(4) 



(5) 



A breakdown of the universality of free fall would manifest 
itself in a particle dependent gravitational acceleration g. 
The solution of Newton's law of motion <j5j reads 



1 



z(t) = z a + v a t - ^gt 2 . 



where zq and vq denote the initial position and velocity of the 
test particle, respectively. Therefore, tests of the universality 
of free fall require identical initial conditions for the two test 
particles. 



So far, we have concentrated on the dynamics of a single par- 
ticle in a linear gravitational field with a well defined initial 
position zq and a well defined initial velocity vq. However, 
in reality it is impossible to prepare the state of the physical 
system with arbitrary accuracy. For this reason we now con- 
sider the dynamics of an ensemble of particles described by a 
classical distribution function 

fo = fo(z,v) . 

The probability to find the particle between z and z + dz with 
a velocity between v and v + dv is given by fo(z, v) dz dv. 
The probability interpretation requires that fo is positive ev- 
erywhere. 

Next we turn to the dynamics of the initial ensemble due 
to a conservative force 

dV 

F(z) = 

Oz 

originating from the potential V = V(z). 

The requirement of conservation of probability leads us 
to the classical Liouville equation 



d d 1 dV d , 

Ot oz m oz Ov ' 



(6) 



subjected to the initial condition fo(z, v) = f(z, v;t = 0). 

For the linear potential Vi, given by ©, this equation 
takes the form 



ot oz ov 1 



(7) 



where we have recalled the definition of specific gravitational 
acceleration ([TJ. 

It is easy to verify, that the solution of (0 reads 

f{z,v;t) = f (z-vt- -gt 2 ,v + gtj . 

This expression brings out most clearly the fact, that the 
dynamics of fo only depends on g, that is on the ratio of grav- 
itational and inertial mass. However, this property does not 
exclude the possibility, that f(z,v;t) can contain in addition 
a dependence on the inertial mass, since the initial distribu- 
tion fo might involve the inertial mass. 

For example, the stationary solution 



f s {z,v) = .yVexp 



U[z) 



2k B T k B T 



(8) 



of the Boltzmann equation representing a gas of colliding par- 
ticles at temperature T in a trapping potential U = U(z) 
involves the inertial mass m. Here, JV and k B denote a nor- 
malization factor and the Boltzmann constant, respectively. 

When we take / s as the initial distribution of our ensem- 
ble of particles propagating in the gravitational field, the final 
distribution f(z,v;t) will obviously involve not only g but 
also the inertial mass m. Hence, in a comparison of the free 



4 



E. Kajari et al. 



fall of two ensembles of particles of different composition, it 
is important to ensure that the initial distributions are identi- 
cal. For the example of two different species of atoms, pre- 
pared in the stationary solution of the Boltzmann equation, 
given by ©, this requirement implies that the temperatures 
and binding potentials have to be adjusted appropriately. 

Hence, in the case of two initial ensembles of particles de- 
scribed by a distribution function in position-velocity space, 
a test of this universality of free fall is more complicated. Al- 
though each member of the ensemble satisfies the universality 
of free fall, the two initial distributions have to be identical in 
order to create two comparable situations. 

3 Wigner function: a few facts 

In classical mechanics kinematics describes motion without 
going into the origin of the motion. On the other hand dynam- 
ics asks for the origin of the motion. In the same spirit quan- 
tum kinematics describes the quantum states and quantum 
dynamics their time evolution. Throughout our paper, this 
distinction will be reflected in the separate treatment of ini- 
tial states and their time evolution. In particular, we consider 
energy eigenstates as natural candidates for initial states. 

To the best of our knowledge the distinction between 
kinematics and dynamics has been spelled out for the first 
time most clearly by Weyl in his book "The Theory of Groups 
and Quantum Mechanics" ll35l . It is interesting that in this 
book, Weyl also defines the concept of averages of symmet- 
rically ordered operators using a distribution function, which 
later became the Wigner function 136U381 . It is this phase 
space function which we use in our quest to analyze how 
the inertial and the gravitational mass manifest themselves 
in quantum mechanics. We devote the present section to a 
brief review of the Wigner distribution and focus on the ele- 
ments most pertinent to the present discussion: the quantum 
Liouville equation and the phase space analog of the time in- 
dependent Schrodinger equation. 

3.1 Definition 

The Wigner function W = W(z,p;t) is a quasi-probability 
distribution which lives in phase space spanned by the posi- 
tion z and its conjugate variable, the momentum p. When the 
state of the quantum system is described by a density operator 
p = p(t), the corresponding Wigner function reads 



and 



W(z,p;t) = 



1 

2nh 



d£e-^/ fi (* + £/2|/}(t)|*-£/2) 



(9) 

where \z) denotes a position eigenstate. 

This expression brings out the fact that the Wigner func- 
tion is real. However, it is not necessarily positive. Moreover, 
the Wigner function satisfies the marginal properties 



dpW(z,p;t) = (z\p(t)\z) = P{z;t) 



(10) 



dz W{z,p;t) = {p\p(t)\p)=P{p;t). 



that is the integrals over the phase space variables p and z 
yield the corresponding quantum mechanical probability den- 
sities P — P(z;t) and P = P(p;t) of the conjugate vari- 
ables. 

The definition (0 of the Wigner function suggests that 
this formulation of quantum mechanics rests on the Schrodinger 
representation and requires a wave function or a density op- 
erator as a starting point. However, this impression is mis- 
leading. The Wigner phase space formulation of quantum 
mechanics is an approach in its own right. In principle, there 
is no need to resort to wave functions or density operators. 
For a more detailed introduction to the Wigner function we 
refer to I 



3.2 Quantum dynamics in phase space 

The dynamics of a quantum state \ip{t)) describing the mo- 
tion of a non-relativistic quantum particle of inertial mass to 
in a potential V = V(z) follows from the Schrodinger equa- 
tion 

d \m) =H\m) (id 



m 



dt 



with the Hamiltonian 



H 



P 

2m 



V(z) 



This description is equivalent to the quantum Liouville 
equation 



d p d dV d 
-7T, + -jr- - tj- tt- - Sf W(z,p;t) 
at to oz oz op 



(12) 



which governs the time evolution of the Wigner function 
W = W(z,p; t). Here the differential operator 



■i 



d 2l+1 V(z) d 2l+1 

Q z 2l+1 Q p 2l+1 



involves only odd derivatives of the potential V and even 
powers of H. 

Although the quantum Liouville equation is in general de- 
rived from the time dependent Schrodinger equation ( fTTT i. we 
could also interpret ( TTZi i as the equation of motion for the 
Wigner function without any reference to the Schrodinger 
formulation of quantum mechanics. Indeed, if we possess a 
priori knowledge about the initial Wigner function, there is 
no need to refer to the time dependent density operator p(t), 
since the Wigner function W(z,p; t) contains all information 
of a quantum system. 
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3.3 Quantum kinematics in phase space 

Energy eigenstates \E) of the time independent Schrodinger 
equation 

H\E)=E\E) (13) 

are the elementary building blocks of quantum mechanics. 
We now briefly motivate the phase space analog of this equa- 
tion. 

Three steps lead to the partial differential equations in 
Wigner phase space for an energy eigenstate: (i) multiply dT3l > 
by (E\, (ii) apply the Weyl- Wigner correspondence IT381 and 
(iii) take the real and imaginary part of the resulting equation. 
In this way we obtain two partial differential equations for the 
Wigner function We = We(z,p) of the energy eigenstate 
that must be mutually satisfied. 

The imaginary part yields the time independent quantum 
Liouville equation 



p d 

m dz 



dVd_ 

dz dp 



je )w E = o 



(14a) 



which through „Sf contains odd derivatives of the potential 
only. 

Equation ( I Hal ) is the quantum Liouville equation ( fT2] i 
with a vanishing time derivative. This feature reflects the fact 
that an energy eigenstate only picks up a phase during its time 
evolution and the Wigner function is bi-linear in the state. As 
a result, this phase factor drops out and the Wigner function 
of an energy eigenstate is time independent. 

From the real part we obtain the analog 



S? )We = EW e (14b) 



h 2 d 2 p 2 
8m dz A 2m 



of the time independent Schrodinger equation which involves 
only even derivatives of the potential in the differential oper- 
ator 

(-1)' (h\ 21 d 2l v(z) d 21 



! 12 



dz 21 dp 21 



It is interesting that both, Jzf c as well as _£f , contain only even 
powers of h. 



where the equal sign holds for pure states. Here we can inter- 
pret the right hand side of this inequality as the effective area 
of phase space taken up by a quantum state 1 39 1 . 

But even if we would choose a normalizable function in 
phase space that satisfies this inequality, it is not clear that 
it represents a physical Wigner function, since it must be re- 
lated to a quantum state with positive semi-definite density 
operator via the Weyl- Wigner correspondence. For a charac- 
terization of the set of all phase space functions that represent 
physical Wigner functions we refer to ll40l . or in the special 
case of Gaussian phase space functions to BTIl . We empha- 
size that the energy eigenvalue equations ( TT4| ) as well as the 
dynamical equation ( fT2b ensure the existence of a valid initial 
Wigner function and its time evolution. 

Since h appears differently in the dynamical and kine- 
matical equations of Wigner phase space, one could imagine, 
at least mathematically, an extended phase space theory of 
quantum mechanics in which both Planck's constants differ 
from each other. For example, a problem where the concept 
of two different h turned out to be useful Il42l is Kramer's 
dilemma and the Langer transformation. However, such an 
extension would probably lead to a physically inconsistent 
theory. In this sense there is a parallelism between the equiv- 
alence of gravitational and inertial mass in general relativity 
and the identity of Planck's constants of dynamics and kine- 
matics in quantum mechanics. 



4 Universality of free fall in Wigner phase space 

Next we consider the partial differential equations ( fT2b and 
(fT4l) determining the Wigner function from phase space for 
the case of a linear gravitational potential given by (0). We 
again analyze quantum dynamics and kinematics separately. 



4.1 Quantum Liouville equation for a linear potential 

For the linear potential © the second and all higher deriva- 
tives vanish. As a result, the operator Jz? , containing all odd 
derivatives of the potential are zero and the equation of mo- 
tion for the Wigner function (fT2l reduces to 



3.4 Constraints on the initial Wigner function 

In order to obtain a unique solution of the quantum Liou- 
ville equation (fT2l . we must specify the initial Wigner func- 
tion Wo = Wo(z,p). However, the choice of Wo is a subtle 
and context dependent enterprise. According to Planck every 
quantum state must take up in phase space at least an area 
2ttH. In the language of Wigner functions this condition as- 
sumes the form 1381 



dz dpW 2 (z,p) 



2nh < 



— oo —oo 



r\ r\ r\ 

+ P- —-mg —)W(z,p;t) = (15) 
at m dz op; 



and no longer involves Planck's constant h. In fact, the quan- 
tum Liouville equation simplifies to the classical Liouville 
equation (Ql for the linear potential when we recall the rela- 
tion 

p = mv (16) 

between the momentum and the velocity of the particle. 

Accordingly, the solution of the quantum Liouville equa- 
tion (15[ is given by 

W(z,mv;t) = W (z-vt-^gt 2 1 m[v+~gt\) , (17) 
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with Wq(z,p) = W(z,p;t — 0) being the initial Wigner 
function of the quantum system. 

The explicit expression ( fTTT i for the time evolution of the 
Wigner function emphasizes again the difference between dy- 
namics and kinematics. Since the quantum Liouville equation 
reduces to the classical Liouville equation in the case of a 
linear potential, each point in phase space propagates accord- 
ing to classical mechanics while maintaining its "weight" as 
given by the initial Wigner function Wo — Wq(z,p). As a 
result, the dynamics only involves the ratio m g /mi. How- 
ever, the initial state may depend on the inertial or gravita- 
tional mass in a nontrivial way and therefore the Wigner func- 
tion W(z,p; t) may not only involve the mass ratio m g /mi. 
Here we find a complete analogy to the classical treatment of 
Sect. 12.21 with one distinct difference: not every initial quasi- 
probability distribution Wa(z,p) in phase space corresponds 
to a possible quantum state. 

We emphasize that the time evolution of the Wigner func- 
tion in a linear potential given by ( TTTb can be related in a 
straightforward manner to the free propagation, as discussed 
in Appendix [A] One important consequence of this consid- 
eration is the fact, that the specific gravitational acceleration 
g does not influence the spreading of the wave packet, but 
only its position along the z-axis. In other words, the vari- 



<- 2 > 



is independent of g. In contrast, the 



expectation value (z) does depend on g. 

Nevertheless, even a classical time evolution can, under 
appropriate conditions, display non-classical features due to 
a non-classical initial state. One prominent example is the 
shrinking [43 44] of a free, radially symmetric wave packet. 



4.2 Wigner function of an energy eigenstate 

Next we turn to quantum kinematics and consider as an ini- 
tial state the energy eigenstate of a quantum particle in the 
linear gravitational potential given by (|4j. The partial differ- 
ential equations determining the corresponding Wigner func- 
tion We = We(z,p) with energy E follow from Eqs. (fl4l . 
In particular, the time independent quantum Liouville equa- 
tion d!4ab reduces to 



p d 

m dz 



d 



(18a) 



whereas the eigenvalue equation ( 1 14bb for the Wigner func- 
tion reads 



dz 2 



8m 



E 



P 

2m 



mgz 



W E (z,p) = 0. 



(18b) 

Since the time independent quantum Liouville equa- 
tion ( 1 1 8ab represents a homogeneous first order partial dif- 
ferential equation, we can apply the method of characteristics 
and deduce that We can depend on the phase space coordi- 
nates z and p only via the classical Hamiltonian 



p 

Hi{z,p) = V mgz 

Am 



(19) 



of a particle in a linear potential. 

The particular functional dependence of the Wigner func- 
tion on Hi(z,p) is then determined by the eigenvalue equa- 
tion ( fT8bl . The full solution 061371 



W E (z,p) = N E ■ Ai 



h 2 mg 2 



H l (z,p)-E 



(20) 



is given in terms of the Airy function Ai = Ai(y), which 
satisfies the ordinary differential equation ||45| 



^-•") Ai| « , = ,) - 



(21) 



Note that the Wigner function We(z,p) is not normalizable. 
As a result the energy E remains a continuous parameter and 
Ne depends on E. 

In Fig.[T]we depict the Wigner function We — We{z,p) 
given by (120) . We recognize a dominant positive-valued ridge 
along the phase space trajectory p = ±p c \(z; E) given by the 
classical momentum 



Pd(z;E) = y/2m(E - mgz) 



(22) 



following from the condition Hi(z.p) = E with the classical 
Hamiltonian ( fT9l . To be precise, at ±p c \ the second derivative 
of the Airy function vanishes. 



W E (z,p) 




Fig. 1 Wigner function We = We(z,p) of an energy eigenstate 
in a linear gravitational potential for E = as given by ( 1201 ). 
The parabolic shape is due to the functional dependence of We on 
the Hamiltonian H9\ . The oscillatory behavior arises from the Airy 
function which follows from the eigenvalue equation d!8bb in phase 
space. 



In the classically forbidden domain of phase space 

mgz 



E < 



P 



2m 



which is inaccessible for a classical point particle moving in 
the linear gravitational potential, the Wigner function decays 
exponentially. 
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In the classically allowed realm of phase space We is os- 
cillatory and can take on negative values. This feature reflects 
the interference nature of quantum mechanics and expresses 
the fact that the energy eigenstate is the superposition of a 
right and a left going wave, as shown in Appendix [C] 

We conclude by noting that the properties of the Airy 
function have also been studied in the context of tunnel- 
ing ll46l . In particular, the dynamics of the Wigner func- 
tion tunneling out of a binding delta function potential in the 
presence of an external static electric field has been studied 

in ma. 



5 Energy wave functions in position space 

The preceding section has employed a phase space analysis 
to reveal the relationship between quantum and classical be- 
havior of a particle in a linear potential. The dynamics of the 
Wigner function are entirely classical, whereas the kinemat- 
ics, and in particular the energy eigenstates, are quantum me- 
chanical. We now address the question how the inertial and 
gravitational mass make their appearances in the correspond- 
ing energy eigenfunctions ue = ue(z) = (z\E) in position 
space. 

In principle, we could obtain the probability density 
u E (z) by taking advantage of the marginal property of the 
Wigner function d20l i together with the integral formula ll48l 

oo 

J d£ Ai(£ 2 +y) = 2^Ai 2 (y/2§). 

— oo 

However, it is equally straight forward to solve the cor- 
responding time independent Schrodinger equation. In the 
present section we pursue this approach and show that the 
spatial modulation of the energy wave function in the linear 
gravitational potential depends on the third root of the prod- 
uct of the gravitational and the inertial mass. 



5.7 Unbounded linear potential 

For a particle in a linear potential, the time independent 
Schrodinger equation reads 



d 2m . 
^2 - [mgz-E] )u E {z) 



0. 



(23) 



Due to the similarity of this equation with the energy eigen- 
value equation in phase space ( I18bb . their solutions must be 
similar in form. Indeed, with the help of the differential equa- 
tion of the Airy function (ISTT i. we can immediately verify that 



u E (z) = -yV E ■ Ai 



mh 2 g 2 



1/3 



[mgz — E] 



(24) 



satisfies d23l . The constant jYe has to be chosen so as to en- 
sure the orthonormality relation 

(E\E') = 8(E — E') 



between two different energy eigenstates \E) and \E') 
It is instructive to cast d24T i into the form 

ue{z) = -JVe ■ Ai (kz — e) , 
by defining the dimensionless energy 



E 



mh 2 g 2 



1/3 



and 



2m 2 g ^ 1/3 

h 2 



(25) 



which has the same physical units as the familiar wave vec- 
tor of a plane wave. For Rb 87 atoms this quantity defines an 
inverse length scale of the order of k w 3.3 x 10 6 m . 

When we insert (Q]) into d25l ). we find that the wave vector 



k = 



2 rrii m g g 

h 2 



1/3 



(26) 



involves the third root of the product of the inertial and the 
gravitational mass. 

Therefore, the spatial modulation of the energy eigenstate 
offers a possibility to compare the masses m g and rrii by a 
method independent of the classical experiments based on 
dynamics. However, it is interesting to note that in the semi- 
classical limit, there is a revival of the universality of free fall 
as shown in AppendixICl 

We conclude by noting that the energy eigenfunction ( 1241 ) 
exhibits a surprising feature when it undergoes a free time 
evolution. In fact, the free propagation of ue does not dis- 
play [49 50 1 the phenomenon of spreading, but just an overall 
acceleration. Recently, such Airy wave packets have received 
great attention in optics and have been realized experimen- 
tally with light 1511521 . For an explanation of this effect in 
Wigner phase space we refer to AppendixlBl 



5.2 Atom trampoline 

Next we insert an infinite repulsive potential wall at z = 0. 
Here, we are not concerned about the nature of this potential 
wall, i.e. whether it originates from electromagnetic forces 
or from gravitational ones that must include the gravitational 
mass. Its purpose is simply to establish a Dirichlet boundary 
condition for the wave function and to provide a Hamilto- 
nian that is bounded from below. Moreover, we emphasize 
that experimental realizations of such a trampoline (or quan- 
tum bouncer) for atoms [ 1 8 - 22 ] , neutrons [ 32 1 , and light ll53ll 
exist. 

The boundary condition 



u E (0) = 



(27) 



on the wave function (l24l i enforces the discrete energy eigen 
values 

' mh 2 g 2 \ 3 



E,, 



ln+1 



(28) 
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Fig. 2 First ten energy eigenfunctions at the corresponding eigen- 
values for the atom trampoline. The potential consists of a linear 
part for < z and a hard wall at z = 0. 



where n = 0, 1, 2, ... and aj denotes the j-th zero [45] of the 
Airy function, as depicted in Fig. 2. 

When we recall the definition (Q} of g the energy eigen- 
values read 

/ 1 \ s 2 _ I 
E n = ( - ftp J mj TOj 3 a n+ i (29) 

2/3 1/3 

and thus depend on m g and m i 1 . 

This is quite a remarkable result because the energy spec- 
trum provides us in principle with a third way to compare the 
inertial and gravitational mass. However, it is not yet clear 
which additional degree of freedom of an atom should be 
coupled to its center-of-mass motion in order to probe the 
energy spectrum with the necessary energy resolution, since 
(m h 2 g 2 /2) 1 / 3 « 2.7 x 1CT 12 eV for Rb 87 atoms. 

6 Conclusion 

In the present paper we have studied the role of the inertial 
and gravitational mass in the quantum mechanical treatment 
of a particle in a linear gravitational potential. Experiments 
involving the dynamics of wave packets, no matter how com- 
plicated the initial state may be, only probe the ratio 



This parameter plays a crucial role in the classical experi- 
ments on the universality of free fall and is closely connected 
to the Eotvos parameter (0. 

However, the probability density of a wave packet might 
well depend on the gravitational and inertial in any arbitrary 
combination. Using the phase space analysis, we can identify 
two sources for this fact, (i) The initial Wigner function is a 
matter of state preparation and thus might involve the inertial 



mass as well as the gravitational mass. For example, when we 
start from a particle in a box, the wave function does not de- 
pend on the mass at all. In contrast, the energy wave function 
of the harmonic oscillator depends on square root of the mass, 
(ii) Due to the marginal property of the Wigner function, the 
time evolved probability density follows from the integration 
of the Wigner function and might create in this way new com- 
binations of the gravitational and the inertial mass. 

A particularly striking example of this additional free- 
dom is provided by the energy wave function of a quantum 
particle in linear gravitational potential. Experiments capa- 
ble of measuring the Airy-function shaped probability density 
would yield information about (mgrrii) 1 ^ 3 . Moreover, spec- 
troscopy of the discrete energy spectrum of a particle bound 
in an atom trampoline would provide us with another scaling 
law (nij/m,) 1 / 3 . 

It is interesting to express these scaling laws in terms of ( 
and the geometrical mean 

M 2 = rrii m g 

of the inertial and gravitational mass. Indeed, we find that 
the spatial modulation of the energy wave function of the lin- 
ear gravitational potential is sensitive to M 2 / 3 , whereas the 
corresponding energy eigenvalues involve the combination 
q1/2 j\f i/3 Tjjg additional information concerning the rela- 
tion between the gravitational and the inertial mass is made 
possible by quantum mechanics and goes beyond the classi- 
cal tests of the universality of free fall. 

In the past century the spectroscopy of the matter wave 
representing the electron in the hydrogen atom has triggered 
the spectacular success of quantum mechanics and quantum 
Electrodynamics (QED). Indeed, the discrete Balmer series 
gave birth to matrix mechanics and the Lamb shift led to 
QED. It would be amusing if in the new century the spec- 
troscopy of matter waves of atoms or Bose-Einstein conden- 
sates would shine some new light on the old question of iner- 
tial and gravitational mass. 
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A Time evolution in a linear potential 

In Sect. 2] we have presented an exact expression for the time 
evolution of the Wigner function in a linear potential. The 
main emphasis of this section was the dependence on the in- 
ertial and the gravitational mass. We dedicate the present ap- 
pendix to represent this dynamics as the product of a shearing 
and a displacement operator acting on the Wigner function in 
phase space. This analysis is reminiscent |54| of the one in a 
time dependent harmonic oscillator. 

We start by recalling ( TT71 ) for the time dependence 



W(z,p;t) = W a [z~^t-ht 2 

* TO 2 



p + mgt 



(30) 



of the Wigner function in terms of the momentum rather than 
the velocity variable. 

In the absence of gravity, that is for g — 0, equation ( f30b 
reduces to 



W f {z,p;t) = W (z-^t,p) 



(31) 



and represents the dynamics of a free particle described by 
the Hamiltonian 

p 2 



The appearance of p in the first argument of the Wigner 
function OTb is responsible for the familiar shearing effect 
of the Wigner function. Due to the marginal property ( [Tol l, 
the shearing in phase space translates into a dispersion of the 
wave packet in position space. 

When we introduce the displacement operator 

9{Z,V)T{z,p) = T{z-Z,p-V) (33) 
and the time dependent shearing operator 

.y(t)T(z 1 p) = T(z-^t,p) (34) 

V TO / 

which act on any phase space function T = F(z,p), we can 
represent the time evolution of the Wigner function in a linear 
gravitational potential ( f30l > in the compact form 

W(z,p;t) = S>(Z l (t),V l (t)),y\t)W (z,p). (35) 

Here we have introduced the time dependent displacement 

1 



(ZiW.Pjft)) = --gt\-mgt 



(36) 



containing the specific gravitational acceleration g. 

We emphasize that the operator 3 provides a representa- 
tion of the Lie group of translations in phase space, whereas 
5? corresponds to the Lie group of shear mappings, accord- 
ingly. In particular, the order of and & is important, as 
reflected by the identity 



y(t)9(Z,V) 



Z + —t,V) y(t). (37) 

TO 



In Fig. [3] we show the time evolution of a Gaussian 
Wigner function Wq = Wo(z,p) in the presence and ab- 
sence of a linear potential confirming the decomposition ( f35T > 
into the product of a time dependent shearing and displace- 
ment. An important consequence of this feature is the fact 
that the linear potential has no influence on the spreading of 
the wave packet. The spreading is solely due to the free time 
evolution of the particle given by the Hamiltonian I 



V(z)k 




Fig. 3 Comparison between the time evolution of a Gaussian 
Wigner function in phase space (bottom) in the presence of a con- 
stant or a linear potential (top), that is for Vf(z) = (blue) or 
Vi (z) = mg z (red). The upper picture also indicates the average 
energy Eq of the wave packet. The center of the Wigner function of 
the free particle propagates along a straight line (blue line), whereas 
the Wigner function in the linear potential follows a parabola (red 
line). The time evolved quasi-probability distributions are depicted 
at three different times and clearly illustrate the fact that the time 
evolution in a linear potential can be decomposed into a free propa- 
gation followed by a shift D = (Zi (t) , Vi (t)) in phase space. 



Another surprising implication of d35l l is the fact that the 
linear gravitational potential cannot influence the interference 
fringes in phase space of a Schrodinger cat state that con- 
sists of a superposition of two Gaussians centered around the 
same position but with slightly different initial momenta. In 
Fig. [4] we show the time evolution of the Wigner function 
corresponding to this superposition state. We note that the in- 
terference fringes do change in the course of time. However, 
this change arises solely from the shearing effect and is not 
due to the gravitational field. Indeed, according to d35l l grav- 
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ity only enters through the displacement of the Schrodinger 
cat obtained from the shearing of free motion. 



Pi 


iHMHti 

► 





z 



Fig. 4 Time evolution in the linear gravitational potential of the 
Wigner function Wo(z,p) corresponding to a Schrodinger cat state 
consisting of a superposition of two Gaussian states located at the 
same position but with different initial momenta. 



The property of the Wigner function to track the classical 
trajectory is not restricted to a time independent linear poten- 
tial, but can be extended to arbitrary g = g{t), resulting in a 
generalized shift (Zi(t), Vi{i)) in phase space. We note that 
this perspective has also been used in the context of many- 
body theory and the Gross-Pi taevskii equation to separate the 
center-of-mass motion from the internal dynamics If55l458ll . 



B Non-spreading wave packets 

An interesting feature of the energy eigenfunction u e given 
by d24l i is the fact [49-52 1 that it does not spread during the 
free time evolution governed by the Hamiltonian H j defined 
by (|32| >. Instead, it preserves its shape and accelerates in pos- 
itive z-direction with a rate g t 2 /2. The representation of the 
quantum dynamics in a linear potential as a product of shear- 
ing and displacement in phase space discussed in AppendixlAl 
offers new insights into the origin of non-spreading Airy-type 
wave packets. 

Indeed, the Wigner phase space description of quantum 
mechanics allows a rather straight forward derivation of this 
effect. We elucidate its connection to a specific invariance 
property of the Hamiltonian dT9b and show that this symme- 
try relation defines a broader class of Wigner functions that 
all correspond to non-spreading wave packets. As an exam- 
ple of this broader class we examine the Wigner function that 
follows from an incoherent superposition of energy eigen- 
states d24l >. We conclude with an alternative view on the effect 



of non-spreading wave packets based on the transformation to 
an accelerated reference frame. 

B.l Equivalence of shearing and displacement 

The Wigner function ( f20b of an energy eigenfunction is time 
independent. Thus, by inserting We into the time evolution 
equation (l35l l. we arrive at the relation 

W B = @3?We- (38) 

When we multiply (l38l ) by the inverse operator of 3l, 
we find 

,yw E = ®- x w E - 

This equation enjoys an interesting interpretation: the shear- 
ing of the Wigner function We is equivalent to a displace- 
ment of We by (—Zi(t), —Vi(t)). Therefore, We preserves 
its shape during the free time evolution given by OTb and 
shifts its position in phase space according to 

W f {z,p;t) = W E {z--gt 2 ,p- mgt} . (39) 

When we take advantage of the marginal property ( fTOb of 
the Wigner function d39b . we find the identity 

oo 

J dpW f (z,p;t) = u E (z-±~gt 2 ), (40) 

— oo 

which indicates that the initial probability density of the en- 
ergy eigenstate (1241) does not spread during the free time evo- 
lution, but accelerates. 

B.2 Broader class of non-spreading wave packets 

The Wigner function We is not the only quasi-probability 
distribution whose marginal P(z;t) exhibits this interesting 
feature. Since the non-spreading behavior of the wave packet 
can be traced back to the invariance relation d38l ). any phase 
space distribution W(z,p) that obeys the identity 

W{z,p) = ®{Z l {t),Vi{t))y{t)W{z : p). (41) 

possesses a probability density P(z; t) in position space that 
does not spread during free time evolution. 

The invariance property ( l4TT i implies that the correspond- 
ing function W(z,p) must depend on the phase space coor- 
dinates z and p only via the classical Hamiltonian 

P 2 

H l {z,p) = ^-+mgz (42) 
2m 

of the particle in a linear potential. 

Indeed, when we apply the operators 3* and 3 to Hi, we 
find with the help of d36i > the relation 

®{Z l {t),Vi{t))y{t)H l {z,p)=H l {z,p). (43) 
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Hence, any phase space function of the form 

W(z,p) = &(H l (z,p)) 



(44) 



automatically satisfies the invariance property (HTt . 

In order to show that the functions d44l i are the only possi- 
ble quasi-probability distributions that satisfy the invariance 
relation fiTT i. we differentiate fiTI ) with respect to the time t 
and insert t = 0. The resulting equation coincides with the 
time independent quantum Liouville equation (I18ab . whose 
general solution is found by the method of characteristics and 
exactly coincides with the class of phase space functions (l44l . 

We emphasize, that the phase space distributions d44b 
cannot be normalized. In fact, when we introduce the new 
phase space coordinates 



3f = Hi{z,p) 



and 



P 1 = P, 



the normalization condition can be expressed as 

oo oo oo oo 

J dz J dpW(z,p) = J dz'^(z') J dp' . (45) 



— oo — oo 



Since the integral over p' diverges, W and therefore its 
marginals cannot be normalized. 

As a consequence, a wave packet that corresponds to a 
Wigner function satisfying the invariance relation (HTl i can- 
not be exactly realized in an experiment and does not allow 
the definition of the expectation values (z) and (p). This fact 
saves the day for Ehrenfest's theorem. 



for the parameters £o = ot{Hi — Eq) and 7 = aa, we obtain 
the exact formula 



x Ai 



a(H l (z,p)-E) + -(aa) i 



(49) 



Here we have introduced the constant 



h 2 mg 2 



The incoherent superposition Wi nc (z,p) of a continu- 
ous distribution g = g(E) of energy eigenstates \E) given 
by ( l46l l differs in general significantly from the Wigner func- 
tion We(z,p) of a single energy eigenstate, as illustrated in 
Fig. [5] In particular, the domains of phase space where the 
Wigner function W lnc assumes negative values have almost 
disappeared. However, due to its functional dependence on 
Hi(z,p), the parabolic profile has survived. 

In the limit a — > for which the energy distribution 
approaches a (5-function centered around Eq, we recover 
from d49l the Wigner function 



lim W inc (z,p) = W Eo (z,p) 

a— ¥0 



of the energy eigenstate |i?o)- 



B.3 Incoherent superposition of energy eigenstates 

Next, we consider an example for the general class of Wigner 
functions ( f44l > with a non-spreading probability density P = 
P(z; t), namely an incoherent superposition 

00 

W inc (z,p)= J dE g(E)W E (z,p) (46) 

— OO 

of Wigner functions We corresponding to a density operator 

00 

p= JdE g(E) \E) (E\ . (47) 

—00 

For the probability distribution g = g(E) to find the par- 
ticle in the energy eigenstate \E), we choose the Gaussian 
distribution 



9(E) = 



1 



V2^ 



: exp 



1 

2^2 



(E - E f 



(48) 



with mean energy Eq and variance a 2 . 

When we recall d20b for the Wigner function We together 
with the integral relation 



7 e -^(^ f 

d£ j=^Ai(0 = e 




Fig. 5 Quasi-probability distribution Wi nc (z,p) corresponding to 
an incoherent superposition of energy eigenstates \E) of a linear 
potential, that is of Wigner functions We(z,p) with a Gaussian 
weight function d48t of energy spread a and average energy Eq. 
In contrast to the Wigner function of a single energy eigenstate, this 

distribution shows strongly suppressed oscillations in the domain 

2 

^— + mgz < Ea and has a rather broad maximum sightly left to 
the classical parabolic phase space trajectory corresponding to Eq. 
It is the average over the Gaussian that has eliminated most of the 
negative contributions in Wigner phase space and has led to a broad- 
ening of the maximum of We- 
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B.4 Transformation to an uniformly accelerated frame 



C Semi-classical limit of an energy wave function 



An alternative explanation for the non-spreading wave packet 
rests on a coordinate transformation [50] from an inertial to 
an uniformly accelerated reference frame. It is worthwhile to 
translate this idea into Wigner phase space in order to obtain 
yet another interpretation of the invariance relation ( |43T >. 

The transformation from the inertial coordinate system to 
the uniformly accelerated reference frame reads 

z' = z - - g t 2 , p' =p- mgt , t' =t, 

with the new phase space coordinates z' and p' and implies 
the following relations on the partial derivatives 



d d d d d d _ d 
dz dz' ' dp dp' ' dt dt' dz' 

As a result, the quantum Liouville equation 

d p d 



dt 



W f (z,p;t) = 



. d 

m9 W 



(50) 



of the free particle expressed in the accelerated reference 
frame takes the form 



m dz dp 1 



dt 



with the transformed Wigner function 

W'{z\p';t') = W f (z' + \~gt! % , p' + m~gt';t') . (52) 

The quantum Liouville equation dSTb is identical to (fT3T > and 
contains a fictitious linear gravitational potential arising from 
the transformation into the accelerated reference frame. 

We now assume that our initial quasi-probability distribu- 
tion for the quantum Liouville equation of a free particle ( f50b 
is given by the Wigner function of an energy eigenstate of a 
linear potential, that is Wf(z,p;0) = We{z,p). The initial 
Wigner function in the accelerated frame governed by dSTb 
follows from d52l for t' = and takes the form 

W'(z',p';0) = W f (z',p';0) = W E (z',p'),. 

Since W E (z',p') satisfies the time independent quantum Li- 
ouville equation ( 1 1 8 at . we conclude from (IBTT l 

^W E (z',p') = 0. 

Hence, the transformed Wigner function exhibits no time evo- 
lution and we find 

W'(z',p';t) = W E (z',p'). 

Insertion of the last expression into d52b finally yields 

W f (z' + ^g t' 2 , p' + mgt'; t') = W E (z',p') , 

which in terms of the original phase space coordinates z and 
p simply reduces to the Wigner function d39l > of the non- 
spreading wave packet. 



The semi-classical limit of quantum mechanics, that is the 
JWKB approximation, has always provided deeper insight 
into the inner workings of quantum theory. In this appendix 
we apply it to demonstrate that the mass ratio m g / rrii emerges 
in the semi-classical limit of the exact energy eigenfunc- 
tion (|24| | far from the classical turning point. Moreover, we 
identify phase space quantization as the origin of the unusual 
scaling law of the energy eigenvalues of the atom trampoline. 



CI Revival of the universality of free fall 
When we recall the asymptotic expansion B31 

Ai(-M) 



1 1 



\y\ 



■cos -\y 



of the Airy function valid for 1 <C \y\, we can approximate 
the energy wave function u E — u E (z) given by (l24l i by a 
superposition 

+A E (z)e~ i ^ E > 



u E (z) = A E {z)e i4 
of two running waves with identical amplitudes 



(53) 



A E {z) 
and opposite phases 



2y / 7r(e — kz)i 



(54) 



Here, p c \ = p c \(z; E) denotes the classical momentum (l22b 
and the turning point z E follows from the condition 



E 



mgz E . 



In this way, the energy wave function u E can be interpreted 
as the most elementary matter wave interferometer. Indeed, 
it consists of a wave running up and one down the linear po- 
tential. Both waves have identical amplitudes and their phase 
difference is governed |59| by an area in phase space deter- 
mined by the strength of the gravitational constant. However, 
the representation ( 1531 is only justified appropriately away 
from the turning point, that is for z < z^. 
From ( l54l > we note the relation 



cX (z:E) = -H 



dz 



As a result we can connect the rate of phase change of the 
JWKB wave given by ( T53b with the classical momentum p c \ 
at the position z. 

When we recall the definition ( 1221 of the classical mo- 
mentum together with the relation ( fTol l between velocity and 
momentum, we arrive at the classical velocity 



= y/2g(zE ~ z) 



(55) 
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of the particle at position z, which only involves g, that is the 
ratio m,g/mi of the gravitational and the inertial mass. Thus, 
in the semi-classical limit we obtain a form of the universality 
of free fall expressed by the relation 



V c \(z) 



dz 



A position-dependent velocity or momentum does not ex- 
ist in the Schrodinger formulation of quantum mechanics. In- 
deed, we either live in position space or in momentum space. 
However, in the semi-classical limit we can obtain mixed 
position-momentum variables. It is interesting that this re- 
sult also follows directly from d20b of the Wigner function, 
as discussed in Sect. 14.21 

However, we emphasize that in the neighborhood of the 
turning point the decomposition ( 1531 does not hold true. As a 
consequence, there the wave function is not able to repro- 
duce (1551 ). where only the ratio of the two masses enters. 
Around the turning point, the energy wave function thus de- 
pends on the third root of the product of the inertial and the 
gravitational mass. 

Our result complements the work of Davies (33). Here, 
the travel time of a wave was calculated based on the Wigner 
time 



dE 



dz 



Pd(z; E) 



dz 



1 



v c i(z;E) 



(56) 



which is valid appropriately away from the turning point of 
the motion. Indeed, the universality of free fall holds true 
again and the travel time only involves the ratio of the two 
masses. 



for the energy eigenvalues. Thus, we find the same prefactor 
as in ( 128) and the approximation ll45l 



for the j-th zero of the Airy function. 

This treatment clearly identifies the origin of the unusual 

\/3 2/3 

scaling laws m i and m 9 of the energy eigenvalues: The 
quantization of the energy levels follows from the quantiza- 
tion of the action, i.e. of an area in phase space. This quantity 
involves both masses. 



D Atomic fountain 

In Appendix [A] we have shown that a constant gravitational 
field has no influence on the phase space interference fringes 
of a Schrodinger cat moving in a linear potential. On the 
other hand atomic fountains provide us with precision mea- 
surements of the gravitational acceleration. In order to bring 
out the similarities and differences between these two mea- 
surement schemes, we first summarize the essential ideas 
that form the basis of an atomic fountain. In the spirit of the 
present paper, which relies almost exclusively on the Wigner 
functions, we then outline the physics of an atomic fountain 
in phase space, thereby sketching the ideas only. In fact, we 
do not intend to develop a complete description of this pre- 
cision instrument. For this purpose we refer e. g. to ll6T1464ll 
and references therein. 



C.2 Phase space quantization as origin of the scaling law 

The form of the eigenvalues in d28b and, in particular, the 
scaling properties of the two masses in ( f29T > can also be de- 
rived in a straightforward manner from the Kramers improved 
Bohr-Sommerfeld rule 

J = j> dz Pcl {z;E) =2irn(n+^\ . (57) 

The Maslov index is 3/4 because there is one "hard" reflec- 
tion at z = contributing 1/2, or equivalently ir to the phase, 
and one soft reflection at the classical turning point contribut- 
ing 1/4, or 7r/2 to the phase 1160) . 

With the help of the classical momentum d22l . we can 
evaluate the action 

4 r^z~ 3 

6 V mg 

which yields, with the quantization condition ((57), the ap- 
proximate expression 




D.l Basic idea 

In the atomic fountain experiments I18I9II , an effective two- 
level atom is moving vertically up against the gravitational 
field of the Earth. Initially, the atom is in its ground state \g). 
However, on its way up a laser pulse prepares a coherent su- 
perposition of its internal levels. The wave vector k of the 
laser is aligned along the z-axis. Since the transition to the 
excited state |e) is accompanied by a momentum transfer due 
to the photon recoil, the atom in the excited state has a differ- 
ent momentum than the atom in the ground state. As a con- 
sequence they accumulate different phases during their prop- 
agation in the gravitational field. After a time r a 7r-pulse 
exchanges the population of the ground and excited state. Fi- 
nally after another time r, a third laser pulse mixes the inter- 
nal levels of the atom and in this way erases the which-way 
information. The quantity measured at the end is the proba- 
bility to find the atom in the ground or excited state. 

The measurement scheme used in the atomic fountain is 
reminiscent of the problem of wave packet interferometry in- 
volving two different molecular surfaces 1165) . Whereas in 
typical molecules these potentials are rather complicated, in 
the atomic fountain they are linear in lowest order. For this 
reason, the latter can be treated fully analytically. 
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D.2 State vector description 



D.3 Determination of the phase difference 



It is straightforward to translate the atomic fountain exper- 
iment into the language of state vectors. We start with the 
initial state 

= W) \g) , 

where \ip) represents the center-of-mass motion along the z- 
axis. After the first laser pulse the atom is in the entangled 
state 



w s T2 



\<P)\g)-ib\i,)\e) 



Here, the unitary operator 



D = e iki 



(58) 



(59) 



imparts the extra momentum from the photon recoil and acts 
as a displacement of the corresponding Wigner function in 
phase space due the interaction between atom and laser. We 
note that we refrain from including additional phases that 
are imprinted on the atom by the laser which must be taken 
into account in a real experiment. In ( T58l . we also assume 
that the interaction between the atom and the laser does not 
otherwise affect the initial center-of-mass wave function \ip), 
which provides a good first order approximation for the real 
situation encountered in atom interferometer experiments. 

The time evolution of the atom in the linear gravita- 
tional field for the time t described by the unitary operator 
JJi = exp(—iHiT/h) yields the state 



1^) = 



1 

71 



The second laser pulse at the time r is a 7r-pulse and inter- 
changes the probability amplitudes of the ground and excited 
state which leads to the expression 



-WUt^)^) - D^UiD^) \g) 



After another period r of unitary evolution, a third laser pulse 
mixes the internal states with a second 7r/2-pulse, and we 
arrive at the final state 



Using the definitions of the momentum displacement opera- 
tor D and the time evolution operator Ui, the two terms in the 
expression for the motional states of the atom \tp g ) or \ip e ) can 
be combined and the nature of the interference of the trajecto- 
ries can be made explicit. With the notation D^ 1 Ui D = U[, 
equation (l60l can be rewritten in the simple form 

As shown below, using the commutation relation between the 
position and momentum operator together with the Baker- 
Campbell-Hausdorff formula, we can prove that 



1 

"71 



Jkgr 



1 



Ui ui |V) , 



which implies the well-known gravity-dependent expression 
for the probability 

P g = Tr{\¥,) <S7| -\g){g\} = ^(l + cos A® 

to find the atom in the ground state. Here, we introduced the 
phase difference A<f> = kgr 2 . 

In order to establish this result, we note that from the ex- 
pansion of Ui, one can prove 

U[ = expi-irt)- 1 Hi D/H) = exp{-iTH[ / h) , 

where the boosted Hamiltonian H[ is just the original Hamil- 
tonian Hi with the momentum p replaced by p + hk, that is 



Hi = D- 1 H,D = 



{p + hk) 2 

2m 



mgz. 



Finally, we use the Baker-Campbell-Hausdorff formula to 
commute U[ and Ui and find 

U{ Ui = Ui U[ e (-"/i) 2 [H!A] = ^ fj, e ikg T 2 j 

which concludes our simple proof of the phase difference be- 
tween the two interfering trajectories. 



where we have introduced the states 



b^UibUi + Uib^Uib 



W (60) 



and 



\^) = -—{%D) 



b- 1 Ui b Ui - Ui b- 1 Ui b 



of the center-of-mass motion of the atom in the ground and 
excited state, respectively. 



DA Phase space description 

In Fig.|6]we represent in phase space the path of the quantum 
state \tp) of the center-of-mass motion in the course of time. 
For an atom that is detected in the ground state \ip g ), we have 
two paths which start from the same point S in phase space 
and end up at the final point E g . One path corresponds to a se- 
quence of unitary evolution, displacement, unitary evolution 
and inverse displacement. Therefore, the atom first moves the 
linear gravitational potential while being in the ground state 
and then after a transition into the excited state continues its 
motion in the potential. It concludes with a transition into 
the ground state at the point E g . The second path starts with 
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the displacement, is followed by unitary evolution and neg- 
ative displacement, and concludes by the inverse displace- 
ment. Here the atom begins its trip by first making a transi- 
tion into the excited state to be followed by motion in the po- 
tential. After de-excitation the atom in the ground state com- 
pletes its path in the gravitational field. Both paths define an 
enclosed area in phase space that determines the phase differ- 
ence between the two amplitudes contributing to \tp g ). 




► Z 



Fig. 6 Phase space representation of an atomic fountain experiment. 
Two 7r/2-pulses separated by the time 2r surround a 7r-pulse at time 
r and cause transitions in the internal states of the atom which are 
accompanied by displacements in phase space along the momentum 
axis. Each parabola section corresponds to the center-of-mass mo- 
tion in the linear gravitational field within the propagation time r. 
The atom starts its journey in the ground state at the point S. Atoms 
eventually measured in the ground state end up at phase space point 
E g , whereas atoms detected in the excited state are to be found at the 
point E e . In both cases two distinct paths lead from the starting point 
S in phase space to the same final point. The area enclosed by the 
two paths expressed in units of h turns out to be twice the phase dif- 
ference A(f> between the two corresponding probability amplitudes. 
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Fig. 7 Atomic fountain experiment represented in position-velocity 
space for typical experimental parameters. Here the momentum 
transfer due to the absorption or emission of a single photon is much 
smaller than the typical momentum associated with the motion in 
the linear gravitational potential. As a result the rather complicated 
closed curve in phase space shown in Fig. [6] simplifies more or less 
to two parabolas {dotted curves) which are connected at z = by 
vertical lines due to the interaction with the laser (solid line in upper 
inset). At the time r, the second laser pulse causes another internal 
transition and a shift in the momentum. Consequently the parabolas 
are interrupted by straight vertical lines representing the momen- 
tum transfer (lower inset). For each trajectory we have marked with 
the symbol r the points in phase space that represent the end of the 
first step of the unitary time evolution. Note that the straight line on 
the inner trajectory (blue dashed line) is therefore traversed thrice, 
whereas the corresponding section on the outer path (purple dashed 
line) is traversed only once. The atom initially in the ground state 
follows the parabola to the point in phase space indicated by r. The 
laser pulse moves the atom vertically up in phase space and the ex- 
cited atom traverses the same vertical line again, but now as part of 
the unitary time evolution in the potential. This picture is an approx- 
imation; in reality the parabolas have a curvature that for parameters 
of the experiment is negligible on the length scale of the momentum 
transfer and invisible even on the magnified scale of the inset. 



Atoms that exit the interferometer in the excited state | ip e ) 
have traversed the same path in phase space as atoms in the 
ground state except that their final point E e on the trajectory 
is different. They either move in the potential, become ex- 
cited, and then follow again the parabola or become initially 
excited, follow the parabola, emit, follow again the linear po- 
tential, and finally get re-excited. 

In Fig.|7]we show the trajectory of the atom in phase space 
during a fountain experiment for realistic parameters ifTTl . 
For this case the vertical straight paths due to the momen- 
tum exchange with the laser pulse are much shorter than the 
parabolas corresponding to the classical motion in a linear 
gravitational potential. 



As a result we can approximate the closed circuit in phase 
space by two parabolas p c i(z; E) and p c \(z\ E + AE) given 
by (|22| |. which are connected by two vertical lines of length 
Ap = hk at z = 0. The difference AE in energy is to first 
order in Ap 

AE = —Ap = — hk , (61) 

m m 

where pe = Pcl(0; E) denotes the momentum at zero poten- 
tial. 

The phase difference Acj) between the two paths turns 
out to be half of the phase space area enclosed by the two 
parabolas (expressed in units of h). In order to provide a sim- 
ple derivation of this result, we first emphasize with the help 
of (T54t that 2h </>(0; E) represents the phase space area inside 
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the parabola p c i (z; E) ranging from z — to the turning point 
Ze- We then denote by 6<fi the area of the circuit expressed in 
units of h and find for small AE the relation 

6<j> = 2 [0(0; E + AE)- 0(0; E)} » 2^AE . (62) 

When we recall the definition (fS6b of the classical transit 
time t from z — to the turning point z^, we find 



= 2-AE. 

h 



(63) 



Finally, substituting Eqs. ( f6TT > into (l63l together with the 
identity = mgr, reduces to 

5<j) = 2fc.gr 2 = 2Z\0 

and proves our conjecture in the limit of small momentum 
transfer Ap < pg. We note that a more thorough analysis 
shows that this result is actually exact. 

This calculation brings to light a remarkable interpreta- 
tion of the atomic fountain experiment. The phase differ- 
ence Atfi between the two paths is half the difference (50 
between the phases of two energy wave functions in a lin- 
ear potential of slightly different energies. Each wave repre- 
sents an interferometer — it is the interference of two counter- 
propagating waves as expressed by (l53l . The phase 20(0; E) 
is the area in phase space enclosed by a straight line at z — 
and the parabola corresponding to the energy E. The area of 
the "crescent moon" is thus the difference of the areas rep- 
resenting the phases of the WKB energy wave functions in a 
linear potential. 



D.5 Comparison between cat and fountain 

We conclude by comparing the interference contained in a 
Schrodinger cat such as the one shown in Fig. |4] and the one 
in the superposition state in the atomic fountain, depicted in 
Fig. [6] In both cases we consider the quantum state of the 
center-of-mass motion in a linear gravitational field. How- 
ever, in the case of the Schrodinger cat state, we face the su- 
perposition of two states that are located at different points in 
phase space. This fact stands out most clearly in the language 
of coherent states \z + ip) of the harmonic oscillator. Indeed, 
the state 



IV'Cat) = ^Cat \z + ip) + \z + ip') 



of the Schrodinger cat corresponds to two Gaussians which 
differ in their momenta p and p' . Here ^Ac&t denotes a nor- 
malization constant. 

In contrast, in the atomic fountain we have the states \ip g ) 
and \ip e ) each of which consists of the superposition of two 
states which only differ in the phase A<\>. In particular, the 
wave function of the center-of-mass motion for the ground 
state reads 



w = ^ 



1 + e 



iAcf> 



where we have introduced the abbreviation ipj = —UiU[\ip). 
The phase turns out to be half the enclosed area in phase space 
shown in Fig. [6] 

Needless to say, the Wigner function of \ip g ) does not 
contain this phase. However, a phase space approach to the 
atomic fountain must contain this phase. It arises when we 
include the internal degrees of freedom into the phase space 
description. Unfortunately, this goes beyond the scope of the 
present paper. 
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